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ABSTRACT 

The detection of outgassing water vapor from Enceladus is one of the great breakthroughs 
of the Cassini mission. The fate of this water once ionized has been widely studied; here we 
investigate the effects of purely neutral-neutral interactions within the Enceladus torus. We find 
that, thanks in part to the polar nature of the water molecule, a cold (~ 180 K) neutral torus 
would undergo rapid viscous heating and spread to the extent of the observed hydroxyl cloud, 
before plasma effects become important. We investigate the physics behind the spreading of 
the torus, paying particular attention to the competition between heating and rotational line 
cooling. A steady-state torus model is constructed, and it is demonstrated that the torus will be 
observable in the millimeter band with the upcoming Herschel satellite. The relative strength of 
rotational lines could be used to distinguish between physical models for the neutral cloud. 

Subject headings: Saturn: magnetospheres — Enceladus — disks 



1. Introduction 

Undoubtedly one of the major discoveries of the 
Cassini mission was the plume of water vapor and 
dust emanating from Saturn's moon Enceladus 



plume emerges (Spitale and Porco 2007 1 possess 



inferred temperatures as high as 180K (Spencer 



(Hansen et al.\2006 


Porco et aZ.|2006 


Waite et al. 


2006 1 . In addition to opening up many lines of re- 



search as to the origin of the vents and their source 



of heat (e.g. 


Roberts and Nimmo|2008 


Meyer and 


Wisdom||2007 


Halevy and Stewart |200 


3 1, the dis- 



covery solved a hitherto outstanding mystery of 
the Saturn system: since 1993 it has been known 
that there exists a large population of hydroxyl 



radicals in the inner magnetosphere (Shemansky 



et al. 1993 Jurac et al. 2002), centered near the 



orbit of Enceladus at 4Rg (IRs = 6.0 x 10 9 cm) 
and extending to R ~ 10-Rg. In order to pop- 
ulate and maintain the OH torus, it was esti- 
mated that Nw ~ 10 28 molecules per second of 
water needed to be injected into the magneto- 



sphere ( Jurac and Richardson 2005 I from a then 



unknown source. Hansen et al. (2006) estimate 



that the release of water from Enceladus is suffi- 
cient (N\y > 5 x 10 27 s _1 ) for the small moon to 
be the source of the observed OH torus. 

The "tiger stripe" regions from which the water 



et al. 2006). Recently, a hotspot temperature 
of 180 K was measured, again by the CIRS in- 
strument (NASA news release 2008-050). The 
gravitational pull of Enceladus (escape velocity 
0.2 km s _1 ) does little to slow the outflow of gas 
at this temperature. Free expansion of water va- 
por at 180K from the surface of Enceladus and 
subsequent intersection of the jet with previously 
emitted vapor will result initially in a narrow 
(AR ~ OARs) torus of vapor with characteris- 
tic temperature < 180K, i.e. velocit y dispersion 
u < 0.5 km s" 1 ( | Johnson ef~aE| 2006 l >- 

The lifetime of the neutral water is limited by 
ionization and dissociation processes whose rates 
are listed in Table [T] |Johnson et al. (2006 1 demon- 
strated that radial spreading of the water torus 
can occur via one of these processes: charge ex- 
change with corotating ions. In this paper we in- 
vestigate a different spreading mechanism. From 
Table[T| it is clear that a neutral molecule does not 
"feel" the influence of the magnetospheric plasma 
or solar UV during a time r > 6 x 10 6 s after pro- 
duction. We shall demonstrate that this is suffi- 



1 



cient time for the neutral water molecules to inter- 
act with each other many times, and for the torus 
to heat and spread to the extent of the observed 
OH torus. 

The dynamics of moon-generated tori have 



been considered by Smyth and Marconi (19931 



in the context of the atomic hydrogen tori of Ti- 
tan and Triton. |Marconi] ( 2003 1 and pecker and| 
Cheng ( 1994 ) investigated the viscous heating of 
the Triton torus via monte carlo simulations. The 
Enceladus torus is unique in consisting of nearly 
pure water vapor (W aite et aZ.||2006| . The partic- 
ular properties of water as a polar molecule make 
the Enceladus water torus a fascinating physical 
and chemical system. Comets are another astro- 
physical system in which almost pure water is 
observed, and much of the work on the spectrum 
and collisions of water vapor has been done in this 



context (e.g. Bensch and Bergm]|2004 Buffa et at. 
2000). We shall make use of these results for the 



cometary case. Global simulations of the water 
torus around Saturn, including all species of neu- 
trals and ions, have been carried o ut by|Jurac et a~L\ 
( [2002] ) and |Jurac and Richardson| fl2005[|2007| ) . In 
these studies, no significant viscous heating was 
found. We discuss the reasons for this discrepancy 
in Section [7] 

In Section |2] we calculate the cross section for 
collisions between water molecules, and estimate 
the collision rate in the water cloud. The nature of 
viscous evolution at low collision rates is discussed 
in Section [3] where we also investigate the impact 
of collisional inelasticity upon both spreading and 
heating rates. Section [4] is devoted to the evalua- 
tion of the inelasticity of neutral water collisions 
in the Enceladus torus, and then in Section [5] we 
find the degree of viscous evolution permitted be- 
fore the neutrals are lost to charge exchange. The 
steady-state column density profile of the heated 
torus is estimated in Section|6] In Section|7]we dis- 
cuss the consequences of the neutral-neutral inter- 
action and make comparisons with previous stud- 
ies. The suitability of the upcoming Herschel mis- 
sion for observing the neutral water is considered 
in Section |8] 

2. Collision cross section 

The water molecule possesses a permanent 
dipole moment /i — 1.85 D — 1.85 x 10~ 18 esu cm, 



Table 1: Rates of external (non neutral-neutral) 
processes. Charge exchange and electron impact 
ionization rates are applicable at radial location 
of Enceladus (4i?s), assuming ions and neutrals 
occupy the same volume. Notes: (a) from Johnson 
et aL| ( |2006| ); (b) from |Burger et ( |2007 i 



Process 




Rate 


(s 


- 1 ) 


Photodissociation 


1.5 


X 


10" 


-7 


(b) 


Photoionization 


7 


X 


io- 


-9 


(b) 


Charge exchange 


1.7 


X 


10" 


-7 


(a) 


Electron impact ionization 


1.6 


X 


10" 


-8 


(b) 


Keplerian orbital 


5 


X 


10" 


-5 




frequency at Enceladus, fix 













which permits long-range interactions between 
molecules and leads to interaction cross sections 
substantially larger than its nominal (~ 5 A ) 
size. We expect the collision cross section a ap- 
propriate for water-water scattering to scale as 
follows. A substantial change of direction occurs 
when the dipole interaction energy is of order the 
molecular kinetic energy at temperature T, i.e. 



3kT 
~2~ 



/'" 



r3 ' 



(1) 



where r is the separation between molecules. We 
then use a ~ irr 2 to obtain 



/2^\ 2/3 



(2) 



t 2 

which for water at 180 K gives a c± 54 A . For 
a more accurate value, appropriate for use in our 
calculations involving viscous spreading, we em- 



ploy the results of Teske et al. (20051, who mea- 



sured experimentally the viscosity of water vapor 
at temperatures 300-450K and a range of densi- 
ties. We use the standard relation for kinematic 
viscosity v, 

u 

(3) 



u 
3na 



where n is the number density and u the thermal 
velocity, to express their results in terms of a colli- 
sion cross section. Their resulting a(T) closely fol- 
lows the power law temperature dependence pre- 
dicted in Eq. [2] Extrapolating to 180K yields 



2 



ct(180) = 80 A 2 . Hereafter we use 



Evolution at low collision rate 



fcoii = 80 



-2/3 



180K 



A 



(4) 



for the water-water scattering cross-section. 

To assess the importance of neutral-neutral in- 
teractions in the Enceladus torus, we perform 
some order of magnitude calculations. Adopti ng a 



production rate of water Nw ~ 1 x 10 8 s (Ju- 



rac and Richardson 2005 1 , and assuming initially 



that the molecules do not interact with each other, 
the total number of water molecules in the torus 
would be 

N w ~ N w t cc (5) 

where r ce ~ 6 x 10 6 s is the lifetime to charge 
exchange with water group ions (see Table [TJ. A 
torus of water at 180K at the orbit of Enceladus 
has extent H ~ 2v/VLk ~ 0.4i?g in the vertical 
and radial directions. We assume that the torus is 
axisymmetric; the azimuthal spreading time ~ 
8ttRs/u ~ 3 x 10 6 s is not dramatically different 
from the lifetime of the individual molecules to 
charge exchange, but we believe axisymmetry to 
be a reasonable approximation for the purposes of 
our investigation. Departures from axisymmetry 
will be considered elsewhere. Dividing the total 
number of molecules by the volume of the torus 
gives a density n~7x 10 4 cm' 3 . We now allow 
ourselves to consider collisions between neutrals. 
At this density, the water-water collision timescale 
is 

Tcoll ~ («cr coll u) _1 ~4xl0 4 S, (6) 

where we used Eq. [4] to obtain cr co u . The neutral 
water molecules can thus collide r cc /r co u ~ 150 
times before suffering a charge exchange collision. 
The rest of this study is devoted to determining 
the consequences of these purely neutral interac- 
tions that occur before the water feels the presence 
of the magnetospheric plasma. Molecular colli- 
sions will lead to viscous spreading and heating 
of the torus, effects which are particularly pro- 
nounced when there is fewer than one collision per 
molecule per orbit (CIkTcoII > !)• For the collision 
timescale derived above, Six t co ii = 2; as the torus 
expands its density will fall and the time between 
collisions will increase, so we expect the torus al- 
ways to be in the low collision rate regime. We 
explore this regime in the next section. 



The action of molecular viscosity in a shearing 
torus is to transport angular momentum outwards 
and to dissipate gravitational energy as heat while 
the center of mass of the torus moves to lower grav- 
itational potential. The viscous behavior of such 
a system depends critically on the number of colli- 
sions per keplerian orbit per molecule, i.e. on the 
parameter £Ikt C o11- When there are many colli- 
sions per orbit, the standard formula for viscosity 
Eq. [3] applies. In this regime a lower collision rate 
results in a higher viscosity, due to the reduction 
in the molecular mean free path. When the col- 
lision rate is low enough, however, (£Ikt co ii > 1) 
the opposite applies. Particles are each "trapped" 
within a radial annulus of width ~ 2eR, where e 
is the orbital eccentricity. In this regime, a de- 
creased collision rate means less contact between 
radial shells and hence lower viscosity. Thus when 
^ifcoii — 1 (as in the Enceladus torus) the vis- 
cosity is in a sense maximal. The viscous heating 
rate per molecule is accordingly also maximal and, 
when there is no cooling, corresponds to an expo- 
nential growth of temperature. This behavior can 
be understood in terms of collisional reorientation 
of the particles' anisotropic velocity ellipsoid in 



the low collision regime (see e.g. Goldreich et al. 
20041. 



The temperature evolution of the torus depends 
critically on whether the viscously dissipated heat 
can be radiated away or otherwise removed from 
translational kinetic energy. Due to their inter- 
nal degrees of freedom, molecular collisions are 
frequently inelastic. To determine the fate of 
the Enceladus torus, we must know the extent 
to which inelasticity affects viscous evolution in 
the low collision rate regime. We quantify the 
problem as follows. Our derivation is based on 



Narayan et al. ( 1994 ) , who work in the shear- 



ing sheet approximation and take moments of the 
Boltzmann equation for the gas particles within 
the sheet. The quantities we wish to study are 
the second velocity moments: the total velocity 
dispersion u 2 and its components along each axis, 
and UrU,),, the shear stress. 



,2. 



Narayan et al. ( 1994 1 assume no spatial gradients 



in these second moments, an approximation only 
strictly valid when u <C W r Q.K, where W r is the 
radial width of our torus. Collisions are taken to 



3 



be isotropic; this is acceptable when we use a col- 
lision cross section derived from measurements of 
viscosity. 



Narayan et al. (19941 considered the steady 



state; we retain time variation (first term in their 
eq. 3.1.11) with an additional parameter £ to rep- 
resent heating; £ is the fractional increase of the 
second moments per collision. The time evolution 
of these moments is then 



du 2 
~dt 



c 

( 

Tool! 



d(u r u^) 
It 



Tcoll 



(u T U^) . (7) 



As in Narayan et al. (19941, the fractional inelas- 
ticity per collision is denoted £. From their Eq. 
3.6.3-3.6.6 and specializing to a Keplerian disk, we 
obtain the following equations for the four non- 
vanishing second moments: 

3(1 + Qu 2 r - \2n K T coXX WH ~ (1 - 0» 3 = W) 

3(i + cK-(i-0« a = W) 

(1 + ~ ^/ffcoii y - 2tt K T coll ul = ((Dl) 

Solving the equations under the assumption 
Qrt ^S> 1 yields 



(12) 

(13) 
(14) 
(15) 

Ux " v ~ UShcr" • (16) 

We see from Eq. [l2]that £ > 0, i- e - the velocity 
dispersion increases with time so long as 



C= (3-110/8, 

u 2 r = ( 20/33) u 2 , 
u\ = (5/33)u 2 , 
u\ = (8/33)w 2 , 
„ (1 + ..a 



lira 



3/11. 



(17) 



Greater inelasticity leads to cooling. In Section [4] 
we determine the inelasticity of water-water colli- 
sions and find that according to the above condi- 
tion, the torus will increase in temperature. 

We can also obtain immediately an expression 
for the viscosity, 



(i-0 



dv<f,/dr 



12 fl 2 K T 



(18) 



where the rate of shear dv^/dr = 3f2i<-/2 in the 
shearing sheet approximation. We see that finite 
inelasticity £ also decreases the viscosity, though 
not as severely as it does the heating rate. Thus 
there can be substantial spreading without heat- 
ing. 

4. Cooling 

4.1. Internal modes of the water molecule 

At the energies considered in this paper, molec- 
ular rotation is the most important internal mode. 
The lowest vibrational mode corresponds to en- 
ergy ~ 0.2 eV. Dissociation of the molecule re- 
quires 5.1 eV, while ionization needs 12.6 eV. 

The water molecule is an asymmetric top with 
principal moments of inertia 1.0, 2.9 and 1.9 x 
10~ 40 g cm 2 . The lowest rotational transitions are 
listed in Table [2] the molecule effectively possesses 
two ground states, corresponding to spin 1 (ortho) 
and spin (para) wavefunctions. 

The population of rotational levels is far from 
thermal equilibrium. Comparing the spontaneous 
emission coefficients from Table |2] to the collision 
time estimated in Eq. [6] reveals At co \\ > 100 (this 
remains true for spontaneous emission from higher 
rotational levels). Thus any rotational modes ex- 
cited by a collision will have decayed radiatively to 
the ground state before the next collision. Energy 
that goes into molecular rotation is effectively lost 
to the system and constitutes cooling. Further, 
when considering collisional excitation we need 
only include transitions from the ground states. 

Other excitation mechanisms 

Internal modes may be excited by other pro- 
cesses, such as solar radiation or microwave back- 
ground pumping, or collisions with electrons (we 



Table 2: The lowest rotational transitions of water, 



from Buffa et al. (2000) 



Transition 

Jk+K- - J' K i +K >_ 


In — 0oo 
(para) 


lio — loi 
(ortho) 


Wavelength (mm) 
Energy (eV) 
A (s- 1 ) 


0.27 
0.005 
5.6 x 10~ 2 


0.54 
0.0025 
3.5 x 10~ 3 



4 



neglect collisions with ions here because we con- 
sider these to be a loss mechanism for the neutral 
water) . Solar pumping results in an excitation rate 



10 s (from values for comets, in Crovisier 



1984 Bcnsch and Bcrgin 2004), still well below the 



spontaneous emission rate. The microwave back- 
ground produces an excitation rate ~ 10~ 7 s _1 . 

The maximum density of electrons in the mag- 
netosphere is n e ~ 100cm -3 (Wahlund et al. 



2005). For any reasonable interaction speed and 



cross section, the rate of electron impacts is much 
less than A > 10 _3 s _1 . We conclude therefore 
that the water molecules are essentially always 
to be found in their ground state, and that ro- 
tational energy cannot be converted into transla- 
tional; rather it can only flow in the other direc- 
tion, via collisions. 

4.2. Transitions between internal modes: 
quantum view 

The typical kinetic energy of a molecule in the 
cold torus at 180K is ^trans = 3fcT/2 ~ 0.023 eV, 
not much larger than the energy spacing between 
the lowest rotational levels. At these lowest tem- 
peratures, we cannot neglect quantum mechanical 
effects. 

The permanent dipole moment of the water 
molecule makes possible long range interactions 
that involve transitions between internal modes. 
In particular, the dipole-dipole interaction po- 
tential leads to the same selection rules as for 



dipole allowed radiative transitions (e.g. Buffa 
et al.: 2000). In the case of rotational transitions, 
this corresponds to A J = ±1,0 (the parity of the 
wavefunction must change, but the total spin must 
not). 

Resonant exchange of a quantum of rotational 
energy is the most common type of interaction 



(Mason and Monchick 19621; since in our case 



molecules are in the ground state when they col- 
lide, these do not concern us. Exact quantum me- 
chanical calculations for the rate of translational 
to rotational energy are unavailable. Instead we 
make a "worst case scenario" estimate of the max- 
imum amount of energy lost per collision, if in ev- 
ery collision a rotation were excited. At 180K, 
the maximum energy loss per collision that still 
permits viscous heating is 0.006 eV (3/11 of the 
thermal energy). 



The ratio of ortho to para water is not known 
a priori. For some comets (temperature ~ 100 K) 
it has been measured at ~ 3 ( Mumma et aZ.|1987] 



Crovisier et al. 1997): we adopt this value. The 



average energy loss to rotation per collision could 
then be at most 0.005 x 0.75 + 0.0025 x 0.25 = 
0.0044 eV, i.e. heating would still occur. In fact, 



Mason and Monchick (1962) estimate that only 



once in every 4 collisions does transfer between 
translational and rotational energy occur, making 
the the average energy lost per collision of order 
0.001 eV. 

4.3. Transitions between internal modes: 
classical view 

From a classical point of view, the water 
molecule consists of frictionless atoms of oxygen 
and hydrogen (two) . Only individual atoms within 
molecules are considered to collide during a colli- 
sion, and they do so instantaneously. This view 
is appropriate when the thermal energy is much 
larger than the energy splitting between states. In 
order to spin up the water molecule, an impactor 
must collide with one of the hydrogen atoms. 
The impactor is an atom within another water 
molecule, with velocity v relative to the molecule 
to be hit. The colliding particle in the impactor is 
most likely to be the (larger) oxygen atom. After 
collision the velocity of the hit hydrogen atom is 



tuqv 
m H + m 



(19) 



If the collider were a hydrogen, vr would be 
halved. The rotation frequency of the impacted 
molecule about the massive oxygen atom is lo = 
vh/L where L is the bond length. The amount of 
energy in rotation is 



2L 2 



(20) 



where Iw is the moment of inertia of the water 
molecule. The ratio of rotational energy to in- 
coming translational energy is given by 



E, 



Iw 
mwL 2 



(21) 



Using the mean moment of inertia I\y = 2 x 
10~ 40 g cm 2 and the bond length L = 1.0 A 2 , we 
obtain 

~ 0.05 (22) 
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for the fraction of energy transferred into molecu- 
lar rotation after a single collision. The fractional 
energy loss to rotation £ = 0.05 is thus substan- 
tially less than the limiting value for heating de- 
rived in Section [3] £u m = 3/11 ~ 0.27. Heating 
can occur. 

We adopt the fractional inelasticity parameter 
derived here £ = 0.05 for the remainder of the 
paper. 

5. Spreading and heating of the torus 

We have established so far that the collision rate 
in the cold torus is conducive to rapid heating, and 
that rotational line cooling is not sufficient to pre- 
vent this heating. As the torus heats and viscously 
spreads, the density will fall and the collision cross 
section will change. To improve upon our sim- 
ple estimates we take these factors into account 
and make a simple model for the time evolution 
of the water torus, assuming it begins unhcated 
and lasts for the charge exchange time, i.e. the 
torus contains iVtot = NT ce , where the source rate 
is believed to be N ~ 1 x 10 28 s _1 . 

We again adopt axisymmetry and assume an 
initial temperature 180K (i.e. u (0) = 0.5 km s _1 ). 
We solve the following simplified differential equa- 
tions. The number density of water molecules in 
the torus centered around 4i?g is 



as a gas with adiabatic index 5/3. We also employ 



n(t) = 



27r(AR s )W r (t)W z (t) 



(23) 



where W r is the full radial width of the torus, and 
W z is the full vertical width. Viscous spreading 
occurs at speed ~ v j AR, where AR is the length- 
scale over which the angular momentum varies. 
Viewing the torus as two segments each of length- 
scale W r /2, the radial spreading rate is approxi- 
mated as 

dw^ 

dt w r {ty K 1 

where the viscosity v is obtained from equation 



(18 1 with £ = 0.05. The velocity dispersion u 



evolves according to 



du 2 

~dt 



u\t) 



311 (*) 



2dn/dt\ 
3 n(t) ) 



(25) 



where the second term in brackets represents cool- 
ing due to expansion of the torus, approximating it 



W z {t) 



2 



8u 2 (t) 



33 



(26) 



for the full vertical width of the torus (twice the 
scale height), where we used Eq. 15 to obtain the 
vertical velocity dispersion. The collision rate is 



n{t)a co \\{t)u(t) , 



(27) 



where a co \\ comes from Eq. [4] Using Eq. [13] for 
the radial velocity dispersion, the torus has initial 
radial width 



W r (0) 



2 J20u 2 (0) 

n^ 4 



33 



(28) 



Figure [T] shows the result of this calculation. 
Although the density and collision rate fall, the 
torus continues to heat. After the charge ex- 
change time, the radial extent of the torus is 
W r ~ 2.5Rg and the vertical scale height is of or- 
der W z ~ 0.7i? 5 ~ 0.18R. Because now W r ~ R 
the neglect of radial variations is no longer valid; 
we address this in the next Section. 




l.xlO" 2.x 10" 3.x 10" 4.x 10° 5.x 10° 6.x 10° 
time (s) 

Fig. 1. — Initial spreading and heating of the 
Enceladus torus, for u(0) = 0.5 km s _1 . Solid 
line: radial width of torus; dashed line: vertical 
scale height. 



We also investigate the effect of varying the ini- 
tial temperature of the water in the torus. At 
lower temperatures, the torus is initially more 



G 




1.x 10" 2.x 10" 3.x 10" 4.x 10° 5.x 10" 6.x 10" 
time (s) 

Fig. 2. — Evolution of radial width of torus, 
with variation in initial temperature. From 
top to bottom, initial velocity dispersion u = 
0.5,0.3,0.1,0.03,0.01 km s _1 . 



dense, and the collision rate may be larger than 
the orbital frequency. In these cases we use the 
following relation for the viscosity 



1 



1 



1 



(29) 



where v\ is the high density-regime viscosity from 
Eq. [3j and comes from the low-collision rate 
expression, Eq. 18 In the high density case, the 



viscous heating rate is small and is set equal to 
du 2 /dt = 2v 1 (3n K /2) 2 . 

Figure [2] displays the results of this calcula- 
tion. When u(0) > Q.l km s _1 the width after 
one charge exchange time is similar to the case 
when u(0) = 0.5 km s _1 . For smaller initial ve- 
locity dispersions, the torus remains in the high 
collision-rate regime and never reaches the critical 
point QkTco\i — 1 where the viscosity is largest. 
The spreading is therefore much less. 

6. Steady-state water torus 

We now adopt the results of the previous sec- 
tion for initial temperature 180K, i.e. 



u z = 0.18 n K R ; u r = 0.28 Q K R, 



(30) 



then calculate the expected steady-state radial dis- 
tribution of water in the torus. In steady state the 



torus, surface mass density S, evolves according to 



as 
at 



3 d 



v{R)Y,R 



1/2 



= — I i? 1 / 2 

RdR\ dRV x ~'~" \) n oss {R) 

(31) 

where the first term on the right is the standard 
equation for radial viscous spreading in a thin Ke- 



plerian disk ( Frank et «L||1992[ ) and the second is 
the loss rate of neutrals to charge exchange and 
ionization. Strictly speaking the above equation 
is valid only when u <C £IkRj i-e. when local disk 
properties can be related to local values of trans- 
port coefficients. Clearly we are close to the edge 
of validity of this approximation. In particular, it 
neglects the high-energy, high-eccentricity tail of 
the particle distribution, and as such we expect 
the amount of spreading calculated to be a lower 
limit. 

From Eq. [18] and [27] the viscosity is given by 
1 u 3 (R)n(R)a con (R) ^ 



12.6 



We convert to surface mass density by means of 

E(fl) = m K20 n(R) ^(R) ' ( 33 ) 

and for the collision cross section we use Eq. [3] 
expressed in terms of velocity dispersion u: 



Ccoll = 80 



0.5 km s 



-4/3 



A 2 



(34) 



Augmenting Eqs. ( 30 1-( 34 ) with Kepler's law 
flu oc i? -3 / 2 and scaling all quantities to the orbit 
of Enceladus, we obtain 



0.29 



R 
Ro 



7/6 



E/m H2 



cm 



cm 2 s 



(35) 



where Rq = 4i?s. 



We then solve the differential equation (31 ) nu- 



merically, setting the surface density equal to zero 
at the edge of the rings, i.e. the rings are con- 
sidered perfect absorbers. We place the edge of 
the rings at an effective radius larger than their 
actual edge at i? c d g c — 2.2Rs- an eccentric orbit 
centered at distance R reaches a minimum radius 
R{\ — e) at pcriapsc, and if the orbit anywhere 
intersects the rings the particle will be absorbed. 
Because we have assigned an eccentricity e ~ 0.28 
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to all particles, the edge of the rings is placed at 
Res = i?cd gc /(l - e) ~ 3.0R S . 

As the torus spreads and heats, the neutral wa- 
ter molecules spend less time in the regions of 
high ion density, and so the lifetime to charge ex- 
change increases. The rate of photodissociation 
to OH+H, however, is independent of radius and 
is similar to the charge exchange rate near Ence- 
ladus (see Table [l]). We consider two cases. First, 
we keep ti oss constant with radius and equal to the 
photodissociation time. This corresponds to treat- 
ing dissociation as a loss mechanism. Second is the 
more physically realistic case, in which we retain 
the OH and treat it in the same way as water: 



OH has a similar dipole moment (see e.g. Sauval 



et al. 1984 1 and mass to water, so its neutral inter- 



actions are similar. We neglect any energy gained 
by the OH upon dissociation, since we expect most 
of this to go to the liberated hydrogen atom. In 
this second case, we adopt the following formula 
for the loss timescale as a function of radius 



'Hoss 



1 + 0.5 



6 x 10 b s, 



(36) 



up to a maximum of 1.8 x 10 s s, corresponding 
to the photoionization timescale. This expression 
was obtained as an approximate fit to fig. 17 of 
Sittler et al. (2008), who calculate the neutral loss 



timescale as a function of radius, including loss to 
charge exchange, electron impact ionization and 
solar photoionization. 

The neutral column density for both cases is 
plotted in Fig. [3] Both profiles are in reason- 
able agreement with the observed OH column den- 
sity from |Jurac et al. \ ([2002]); |Jurac and Richard 



son 



(2005) 



The true OH column is bracketed by 
the two lines plotted. It must be noted that a 
direct comparison between these curves and the 
measured OH column is inappropriate for the fol- 
lowing reasons: first, we have plotted water or 
water+OH, not OH alone; and second, we have 
neglected spreading of the torus via the other es- 



tablished mechanism, charge exchange (Johnson 

It is of interest to compare the derived full 
width of the torus with the result of the calculation 
in Section [5] From Fig. [3] the full width between 
points with Nw = Mv,max/e is AR ~ 2.7 Rs, con- 
sistent with W r ~ 2.5Rs derived in Section [5] 



> 12.5- 




Fig. 3. — Calculated column density profile of the 
Enceladus water torus. Dashed line: treating pho- 
todissociation as a loss mechanism. Solid line: 
treating only ionization and charge exchange as 
loss mechanisms; line represents water + OH col- 
umn density in this case. 



7. Implications and comparison with pre- 
vious work 

Our results serve to show that significant 
spreading of the torus can occur via neutral- 
neutral interactions alone, and thus such inter- 
actions, with correct cross sections, should be in 
future included in the global models such as those 



of Jurac and Richardson (2005) which include all 



relevant processes. 

A further issue that can be addressed is depo- 
sition of water on to the outer parts of the rings. 



Jurac and Richardson ( 2007 1 estimate using their 



model that 17% of the water from Enceladus is 
precipitated on to the rings. Farrell et al. ( 2008 1 



report measurements of a decrease in plasma den- 
sity at the edge of the rings, and use this to cal- 
culate the rate of absorption of ions. In both of 
the cases shown in Fig. [3] about 70 kg s _1 (25%) 
of the total 280 kg s _1 released by Enceladus pre- 
cipitates on to the rings. This is probably a lower 
limit, both because other processes also contribute 
to loss on to the rings; and because our model 
probably represents a lower limit to the degree of 
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spreading of the torus. 

The detailed models of |Jurac et al. (20021; Ju- 



rac and Richardson] ([2005 2007) do not show 



significant viscous heating or interactions be- 
tween neutrals. We suggest that this is because 
their adopted water-water cross section was that 



of atomic oxygen (Jurac and Richardson 20051, 
which has radius r ~ 0.6 A. In fact, the character- 
istic collisional size of the molecule is substantially 
larger, at least that of the 2 ~ ion (r ~ 1.4 A), 
and is further enhanced by the permanent dipolc 
moment as detailed in the Introduction. Because 
there are fewer collisions between neutrals before 
they are impacted by other processes, they do not 
spread or heat significantly in the models of Jurac 
et al. 



Marconi (2003) performed monte carlo simu- 



lations of the low-collision rate atomic hydrogen 
torus around Neptune's moon Triton. This ap- 
proach follows individual particles and so does not 
suffer the potential inaccuracies of our thin disk 
model. The comparison between Marconi's results 
and ours, however, is encouraging. After a similar 
number of collision times, the torus has spread by 
only a slightly larger amount than we find in our 
Section [6j 

8. Observing the water torus remotely 

The low density of the extended water torus 
makes it difficult to observe in situ and remotely. 
Only in the fluorescent emission of the hydroxyl 
radicals has the extended neutral torus been seen 



(Shcmansky et al. 1993 Jurac et al. 20021, while 



only very close to Enceladus are the neutral den- 



sities large enough to be measured in situ (Waite 
et oZ. 1120061) . Once ionized, detection is easier, but 



the additional electromagnetic forces felt by ions 
and electrons renders the inference of the distri- 
bution of their parent neutral population difficult, 



although this has been attempted (Sittler et al. 



2008 ). The majority of the extended neutral water 
cloud has never been directly detected. The up- 
coming Herschel mission has the ability to change 
this. 

One of the science goals of the Herschel Space 
Observatory, due for launch in 2008, is the obser- 
vation of water rotational lines in the nearby uni- 



verse (e.g. de Graauw et al. 20061 with unparal- 
leled sensitivity and angular resolution. The HIFI 



instrument on board Herschel provides continuous 
frequency coverage from 480 to 1250 GHz, cov- 
ering ~ 10 of the lower rotational transitions of 
the water mol ecule, at frequency resolu tions up to 
vjAv = 10 7 ( |de Graauw et a?.||2006) . The 3.5 
m diameter of Herschel provides angular resolu- 
tion at 560 GHz of ~ 0.5 arcmin, corresponding 
to ~ 3.6i?5 at the distance of Saturn. This resolu- 
tion is sufficient to crudely map the neutral torus. 

We now estimate the signal strength from the 
torus. For definiteness, we consider the Iiq — loi 
(ortho) transition at 560 GHz, and use our steady- 
state hot torus model from the previous section. 
We take the optimistic view that every collision 
leads to an excitation of this line (but neglect emis- 
sion due to electron collisional and solar excita- 
tion, which will be significant), and calculate the 
flux per beam of Herschel as a function of distance 
along the equator, where Saturn is assumed to be 
at equinox. 

The received flux is expressed in the standard 
way in terms of an antenna temperature T ant : 



\ 2 B V 



(37) 



where B v is the specific intensity of the source 
(e.g. Burke and Graham-Smith 20011. Emission 
in the 560 GHz line is spread by both thermal and 
orbital Doppler shifts; we adopt a linewidth char- 
acteristic of the local velocity dispersion, which 
is a fraction of the orbital velocity in our model, 
Av = Udi Bp (i£)/c. The typical linewidth is then of 
order Av ~ 20 MHz, i.e. v/Av ~ 3 x 10 4 . 

Integrating along lines of sight through the 
equatorial plane of the torus, and neglecting beam 
dilution due to vertical variation in density, we ob- 
tain the antenna temperature profile shown in Fig. 
[2] With a system temperature T sys = 160 K in the 
lowest band of the HIFI instrument, Herschel has a 
5-sigma detection limit of 1.6 mK, in 1 hour of ob- 
servation at frequency resolution 10 4 (de Graauw 
et al. 2006 1, and this is also plotted in Fig. |4| 



MS 



the horizontal dotted line. 

Under the assumed conditions, then, the neu- 
tral water torus is readily observable by Herschel. 
The angular resolution prevents the detailed ex- 
amination of the spatial distribution, but by com- 
paring the flux in various lines it should be possible 
to determine the collision rate in the torus and to 
gauge the temperature of the water via linewidths. 
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Fig. 4. — Antenna temperature in the lowest 
ortho- water rotational line, with beam centered 
on equator. Horizontal axis shows distance along 
equator from center of planet. Signal integrated 
through cloud assuming it is optically thin. Hor- 
izontal dotted line shows 5-sigma detection limit 
for Herschel HIFI instrument after 1 hour of inte- 
gration at frequency resolution vjAv = 10 4 . 



A more comprehensive model of the torus should 
be used to make detailed predictions for different 
physical models. 

8.1. Absorption by the torus 

In the above the torus was assumed optically 
thin in the 560 GHz line. In fact the optical depth 
is close to unity and a bright background source 
could allow observation of the torus to further dis- 
tances from the planet, where less emission occurs. 

The optical depth in the line is given by (e.g. 
Rybicki and Lightmanfl979 1 



Tin 



Nw^ 2 ^4iio-i 
Av 



(38) 



where Av is the linewidth and N\y is the line 
of sight column density in water molecules in the 
lower of the two rotational states (the vast major- 
ity of molecules, as discussed in Section |4| . Using 
the spontaneous emission coefficient from Table [2] 
and the line of sight column density at 6 Rs from 
our steady-state model in Section [6] we obtain 



The torus is thus readily detectable in absorption 
against any sufficiently bright millimeter source 
that passes behind it. It small angular size, how- 
ever, makes the existence of any such source rather 
unlikely; the cosmic microwave background could 
be used, but because of its near isotropy on the 
sky, the strength of the signal would be reduced 
by a factor ~ vn sp /c. 

It is clear from this work that neutral-neutral 
interactions in the Enceladus torus cannot be ne- 
glected; they alone can account for the observed 
radial extent of the hydroxyl torus, though un- 
doubtedly spreading also occurs through charge 
exchange reactions. The water- water collisions, 
with the cross sections calculated here, should be 
included in future global simulations of the torus. 
This may impact our estimates of ice deposition 
rates on the rings, as well as the required source 
rate to sustain the observed torus. The Herschel 
mission will provide a timely means of observing 
the neutral water for the first time, by means of 
its rotational lines, and may be able to distinguish 
between alternative neutral cloud models. 

The author wishes to thank R. Narayan and P. 
Goldreich for productive conversations. 
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